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Abstract 

We provide evidence in favor of the conjectured duality between color and kinematics for the 
case of nonsupersymmetric pure Yang-Mills amplitudes by constructing a form of the one-loop 
four-point amplitude of this theory that makes the duality manifest. Our construction is valid 
in any dimension. We also describe a duality-satisfying representation for the two-loop four- 
point amplitude with identical four-dimensional external helicities. We use these results to obtain 
corresponding gravity integrands for a theory containing a graviton, dilaton, and antisymmetric 
tensor, simply by replacing color factors with specified diagram numerators. Using this, we give 
explicit forms of ultraviolet divergences at one loop in four, six, and eight dimensions, and at two 
loops in four dimensions. 
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I. INTRODUCTION 



Recent years have seen remarkable progress in computing and understanding scattering 
processes in gauge and gravity theorieSjboth for phenomenological and theoretical applica- 
tions. (For various reviews see refs. [j], Q.) In particular, various new structures have been 
uncovered in the amplitudes of these theories (see, for example, ref. (3|). One such structure 
is the duality between color and kinematics found by Carrasco, Johansson, and one of the 
authors (BCJ) 0, @]. This duality is conjectured to hold at all loop orders in Yang-Mills 
theory and its supersymmetric counterparts. Besides imposing strong constraints on gauge- 
theory amplitudes, whenever a form of a gauge-theory loop integrand is obtained where the 
duality is manifest, we obtain corresponding gravity integrands simply by replacing color 
factors by specified gauge-theory kinematic numerator factors. 

The duality between color and kinematics has been confirmed in numerous tree-level 
studies 0-11], including the construction of explicit representations for an arbitrary number 
of external legs [l2j]. At loop level, the duality remains a conjecture, but there is already sig- 
nificant nontrivial evidence in its favor for supersymmetric theories [B, 13-3 and for special 
helicity configurations in nonsupersymmetric pure Yang- Mills theory 0, 18| . Here we pro- 
vide further evidence in favor of the duality at loop level, explicitly showing that it holds for 
pure Yang-Mills one-loop four-point amplitudes for all polarization states in D dimensions. 
We also show that a previously presented form of the two-loop four-point identical-helicity 
amplitude of pure Yang-Mills theory [l9| satisfies the duality without rearrangements once 
an additional diagram that integrates to zero is added. Ref. (Ef already noted the existence 
of a form that satisfies the duality; the form we use here is slightly different with nonlocal 
numerators. 

In order to construct the one- loop four-point pur e Yang-Mills amplitude, we use a D- 
dimensional variant [20| of the unitarity method j21] . Our construction begins by finding 
an ansatz for the amplitude constrained to satisfy the duality. Since the amplitude is fully 
determined from its D- dimensional unitarity cuts, we obtain a form of the amplitude with the 
duality manifest by enforcing that the ansatz has the correct unitarity cuts. The existence 
of such a form where both the duality and the cuts are simultaneously satisfied is rather 
nontrivial. We do not use helicity states tied to specific dimensions but instead use formal 
polarization vectors because we wish to have an expression for the amplitude valid in any 
dimension and for all states. The price for this generality is that the expressions are lengthier. 
Since the constructed integrand has manifest BCJ duality, the double-copy construction 
immediately gives the corresponding gravity amplitude in a theory with a graviton, dilaton, 
and antisymmetric tensor. 

We use these results to study the ultraviolet divergences of the corresponding gravity 
amplitudes. Recent years have seen a renaissance in the study of ultraviolet divergences in 
gravity theories, in a large measure due to the greatly improved ability to carry out explicit 
multiloop computations in gravity theories 0, [l3-15, 22h24]. The unitarity method also 
has revealed hints that multiloop supergravity theories may be better behaved in the ul- 



traviolet than power-counting arguments based on standard symmetries suggest [25|. Even 



pure Einstein gravity at one loop exhibits surprising cancellations as the number of external 



legs increases [26|. The question of whether it is possible to construct a finite supergravity 
is still an open one, though there has been enormous progress on this question in recent 
years, including new computations and a much better understanding of the consequences 
of supersymmetry and duality symmetry (see e.g. refs. 27|, |28j). In half-maximal super- 
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gravity 2t|, two- and three- loop examples are known where the divergence vanishes, yet the 
understanding of the possible symmetry behind this vanishing is incomplete [23|, 24, 28, 30 . 
The duality between color and kinematics and its associated double-copy formula offer a 
new angle on the ultraviolet divergences in supergravity theories jE], [l3|, 0, 24, 31]. Here 
we explore the ultraviolet properties of nonsupersymmetric gravity from the double-copy 
perspective. 

We use the gravity integrands constructed via the double-copy property to determine the 
exact form of the ultraviolet divergences. We do so at one loop in dimensions D = 4,6, 8. The 
ultraviolet properties of one-loop four-point gravity amplitudes have already been studied in 
some detail over the years, including cases with scalars or antisymmetric tensors coupling to 



gravity |26l. I32H36I]. so no surprises should be expected, at least at four points. Nevertheless, 
it is useful to look in some detail at the ultraviolet properties to understand them from the 
double-copy perspective. Here we examine the four-point amplitudes in a theory of gravity 
coupled to a dilaton and an antisymmetric tensor, corresponding to the double copy of pure 
Yang-Mills theory. While related calculations have been carried out, we are unaware of 
any calculations of the ultraviolet properties in the theory corresponding to the double-copy 
theory. 

We find that in D = 4, there are no one-loop divergences in amplitudes involving external 
gravitons, though there are divergences in the remaining amplitudes involving only external 



dilatons or antisymmetric tensors, as expected from simple counterterm arguments [32 
By two loops, even the four-graviton amplitudes contain divergences, as we demonstrate 
by computing the form and numerical coefficient of the divergence. In the two- loop case, 
the divergence is proportional to an R 3 operator which gives a divergence in the identical- 
helicity four-point amplitude. In D = 6 and D = 8, we find one-loop divergences in the 
four-external-graviton amplitudes. These results are not surprising and are in line with the 
earlier studies. Our conclusion is that, by itself, the double-copy structure is insufficient 
to render a gravity theory finite in D = 4 and requires additional ultraviolet cancellations, 
such as those from supersymmetry. 

This paper is organized as follows. In Section [Til we briefly review the duality between 
color and kinematics and the double-copy construction of gravity. In Section IIII[ we present 
the construction of the duality-satisfying pure Yang-Mills numerators at one and two loops. 
Then in Section HV] we study the ultraviolet properties of gravity coupled to a dilaton and an 
antisymmetric tensor at one loop in four, six, and eight dimensions. In the same section, we 
also present the ultraviolet properties at two loops in four dimensions. Finally, in Section [V] 
we give our conclusions. Appendices evaluating two- loop integrals needed in Section lTV Bl are 
included. Appendix [A] focuses on extracting the divergences in dimensional regularization. 
This procedure mixes infrared and ultraviolet divergences; so, in Appendix [B] we give the 
infrared divergences that must be subtracted to obtain the ultraviolet ones. Appendix [C] 
evaluates the integrals using an alternative method for obtaining the ultraviolet divergences 
more directly, by introducing a mass to separate out the infrared divergences from the 
ultraviolet ones. 
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FIG. 1: The basic Jacobi relation for either color or numerator factors. These three diagrams can 
be embedded in a larger diagram, including loops. 

II. REVIEW OF BCJ DUALITY 

An m-point L-loop gauge-theory amplitude in D dimensions, with all particles in the 
adjoint representation, may be written as 

jL-\oo P _ ,L m-2+2L V V / TT ^ 1 ^ (2 1) 

where g is the gauge coupling constant. The first sum runs over the m\ permutations of the 
external legs, denoted by S m . The Sj symmetry factor removes any overcounting from these 
permutations and also from any internal automorphism symmetries of graph j. The j'-sum is 
over the set of distinct, nonisomorphic, m-point L-loop graphs with only cubic (i.e., trivalent) 
vertices. These graphs are sufficient because any diagram with quartic or higher vertices can 
be converted to a diagram with only cubic vertices by multiplying and dividing by the appro- 
priate propagators. The propagators appearing in the graph are 1 / Yl a .p&.. The nontrivial 
kinematic information is contained in the numerators rij and depends on momenta, polariza- 
tions, and spinors. In supersymmetric cases it will depend also on Grassmann parameters, 
if a superspace form is used. The loop integral is over L independent D-dimensional loop 
momenta, p\. Finally, Cj denotes the color factor, obtained by dressing every vertex in graph 
j with the group-theory structure constant, f abc = iy2f abc = Tr( [T a , T b ]T c ), where the 
hermitian generators of the gauge group are normalized via Ti(T a T b ) = 5 ab . 

The numerators appearing in eq. (12.11) are by no means unique because of freedom in 
moving terms between different diagrams. Utilizing this freedom, the BCJ conjecture is 
that to all loop orders, representations of the amplitude exist where kinematic numerators 
obey the same algebraic relations that the color factors obey [3, H|. In ordinary gauge 
theories, this is simply the Jacobi identity, 

Ci = cj - c k rii = rij - n k , (2.2) 

where i, j, and k label three diagrams whose color factors obey the Jacobi identity. The 
basic Jacobi identity is displayed in Fig. [TJ The identity generalizes to any loop order with 
any number of external legs by embedding it in larger diagrams, where the other parts of 
the diagrams are identical for the three diagrams. Furthermore, if the color factor of a 
diagram is antisymmetric under a swap of legs, we require that the numerator obey the 
same antisymmetry, 

Ci ->■ -Ci =>- n» ->> -rii . (2.3) 
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The duality was noticed long ago for four-point Feynman diagrams 37j; beyond this, it is 
rather nontrivial and no longer holds for ordinary Feynman diagrams. We note that the 
numerator relations are nontrivial functional relations because they depend on momenta, 
polarizations, and spinors, as discussed in some detail in refs. 0, [U 13]. 

While a complete understanding of the duality and its consequences is still lacking, a 
variety of studies have elucidated it, especially at tree level. In particular, this dualit 
leads to nontrivial relations between gauge-theory color-ordered partial tree amplitudes 
The duality (12. 3p has also been studied in string theory p, 0, 39]. In the self-dual 

. Explicit 
Although 
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representation 40 



dimensional Lie algebra 



case, light-cone gauge Feynman rules have been shown to exhibit the dualit; 
forms of n-point tree amplitudes satisfying the duality have been found |12| 
we do not yet have a complete Lagrangian understanding, some progress in this direction 
can be found in refs. (H nj. The duality (12. 2p does not need to be expressed in terms 
of group structure constants but can alternatively be expressed in terms of a trace-based 
Progress has also been made in understanding the underlying infinite- 
, |4l| responsible for the duality. The duality between color and 
kinematics also appears to hold in three-dimensional theories based on three algebras [42], 
as well as in some cases with higher- dimension operators j43[ . Some initial studies of duality 
and its implications for gravity in the high-energy limit have also been carried out (44| . 

At loop level, the duality remains a conjecture, but there is already nontrivial evidence 
in its favor, especially for supersymmetric theories. At present, the list of loop- level cases 
where duality satisfying forms of the amplitude are known to hold includes: 

• Up to four loops for four-point M = 4 super- Yang- Mills 0, [l3| in a form valid in D 
dimensions; 



up to two loops for five external gluons in M = 4 super- Yang-Mills theory (45| ; 
up to seven points for one-loop amplitudes in M = 4 super- Yang-Mills theory (46| ; 
up to two loops for four-point identical- helicity pure Yang-Mills amplitudes [EJ; 
throu gh n points for one-loop all-plus- or single-minus-helicity pure Yang-Mills ampli- 



tudes 18 



through four loops for a two-point (Sudakov) form factor in = 4 super- Yang-Mills 
theory 0; 



• one-loop four-point am plit udes in Yang-Mills theories with less than maximally super- 
symmetric amplitudes [171 ]. 

In this paper, we add the nonsupersymmetric pure Yang-Mills one-loop four-point ampli- 
tude in D dimensions to this list. Besides direct constructions, we note that the duality 
also appears to be consistent with loop-level infrared properties of both gauge and gravity 
theories [l6j . 

The most significant consequence of the duality is the ease with which gravity loop 
integrands can be obtained from gauge-theory ones, once the duality is made manifest 0, |5( . 
One simply replaces the color factor with a kinematic numerator from a second gauge theory, 



(2.4) 
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FIG. 2: The Jacobi relations determining either color or kinematic numerators of the four-point 
diagrams containing either a triangle or internal bubble. 



This immediately gives the double-copy form of gravity amplitudes, 
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(2.5) 



where rij and rij are gauge-theory numerator factors. An important aspect of the double- 
copy formula is that only one of the two sets of numerators needs to satisfy the duality 
The double- copy formalism has been studied at loop level in some detail in a 
variety of cases 0, 0, M, EH M, 0, H S . 



III. CONSTRUCTION OF DUALITY-SATISFYING INTEGRANDS 

We now describe the construction of a duality-satisfying representation of the one-loop 
four-point amplitude in pure Yang-Mills. Since we want the form to be valid in all dimen- 
sions and for all D — 2 gluon states, we use formal polarizations instead of helicity states. 
This complicates the expression for the amplitude, but has the advantage that it allows us 
to straightforwardly study the amplitude and its gravity double copy in various dimensions. 
In this section, we also show that the form of the two-loop pure Yang-Mills identical-helicity 
amplitude given in ref. [19| satisfies BCJ duality once a contribution that integrates to zero 
is added. In Section HVBt we use this amplitude to show that although four-graviton am- 
plitudes are ultraviolet finite in D = 4 at one loop, they diverge at two loops, in accordance 
with expectations. 



A. One Loop 

For a one-loop n-point amplitude, the duality (12. 2 p can be used to express kinematic 

numerators of any diagram directly in terms of n-gon numerators. In particular, for the 
four-point case we have two basic relations determining triangle and bubble contributions 
from box numerators as illustrated in Fig. [2], 

n 12(34);p — ^1234;p — "^1243;p 5 
™(12)(34);p = ^12(34);p — ™21(34);p • (3.1) 
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FIG. 3: The color or kinematic Jacobi relations involving a bubble on an external leg or a tadpole. 
These diagrams have vanishing contribution to the integrated amplitude. 

The labels 1,2,3,4 refer to the momenta and states of each external leg, while the label p 
denotes the loop momentum of the leg indicated in Fig. [2j (The parentheses in the subscript 
of the numerators indicate which external legs are pinched off to form a tree attached to 
the loop.) Note that in the figure the momentum of each internal leg of each diagram is 
the same as in the other two diagrams except for the single internal leg that differs between 
the diagrams. In general, the bubble and triangle contributions are nonvanishing; indeed, 
this explicitly holds for the BCJ representation of the one-loop four-point amplitude of pure 
Yang-Mills theory that we construct. 

Besides the diagrams in Fig. [2], there are diagrams with a bubble on an external leg and 
diagrams with a tadpole, as shown in Fig. [3j The duality also determines the numerators of 
these diagrams via 

n l(234);p = ^12(34);p + ^l(43)2;p , 
n (12U);p = ^(12)(34);p ~ ™(12)(34);-p , 

n (i234);p = ™l(234);p ~ ™l(234);-p , (3.2) 

corresponding respectively to the three relations in Fig. [3J (On the final line in eq. (13.21) . 
the hat marks leg 1 as the location where the tadpole is attached.) We use these equations 
to impose the auxiliary constraint that the tadpole numerators determined by BCJ duality 
vanish identically and that all terms in the bubble-on-external-leg diagrams integrate to zero 
as they do for Feynman diagrams. Thus, these diagrams are not necessary for determining 
the integrated amplitudes (though in D = 4 the bubble-on-external-leg diagrams do affect 
the Yang-Mills ultraviolet divergence). 

Once we impose the BCJ conditions, the amplitude is entirely specified by the box nu- 
merators. Our task is then to find an expression for the box numerators such that we obtain 
the correct amplitude. It is useful to impose a few auxiliary constraints to help simplify the 
construction: 

1. The box diagrams should have no more than four powers of loop momenta in the 
pure Yang-Mills case, matching the usual power count of Feynman-gauge Feynman 
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FIG. 4: The two-particle cut with a massless invariant. An off-shell continuation of such cuts is 
generically required. The exposed intermediate legs are on shell. 



diagrams. 



2. Each numerator written in terms of formal polarization vectors respects the symmetries 
of the diagrams. In particular, this condition implies that once a box diagram with 
one ordering of external legs is specified, the other orderings are obtained simply by 
relabeling. 

3. The numerators of tadpole diagrams vanish prior to integration. 

4. All terms in the bubble-on-external-leg diagrams integrate to zero, in the same way 
as they do for Feynman diagrams. 

While it is not necessary to impose these conditions, they greatly simplify the construction. 
They ensure that the type of terms that appear in the ansatz are similar to those of ordi- 
nary Feynman-gauge Feynman diagrams, avoiding unnecessarily complicated terms. (Using 
generalized gauge invariance, one can always introduce arbitrarily complicated terms into 
amplitudes, which cancel at the end.) 

The first three conditions simplify the construction by restricting the number of terms 
that appear. The purpose of the fourth auxiliary constraint is a bit more subtle. This 
constraint eliminates terms which can complicate the extraction of ultraviolet divergences 
in both Yang- Mills and gravity. In Feynman gauge, the bubble-on-external-leg contributions 
have the property that they are proportional to (k 2 )( D ~ 2 ^ 2 . In D > 4, after canceling the 
intermediate on-shell propagator, the diagrams are proportional to {k 2 )^ D ~ 4 ^ 2 = and hence 
vanish. Near D = 4 the ultraviolet divergences mix with infrared divergences which need 
to be subtracted. In gravity, there are extra powers of numerator momenta which ensure 
that even in D = 4 there are no contributions arising from bubbles on external legs. The 
fourth condition ensures that we can treat the bubble-on-external-leg contributions in the 
same way as for Feynman diagrams. If we do not impose this condition, we can obtain 
additional nonvanishing contributions that are problematic to analyze because they contain 
an intermediate on-shell propagator, which would need to be regularized in some way by an 
off-shell continuation. While these can, in principle, be analyzed by considering a regularized 
unitarity cut of the type shown in Fig. HJ it is best to avoid this complication. 

We start the construction with an ansatz containing all possible products of £j • Ej, p ■ £j, 
ki ■ Ej, p ■ ki, p ■ p, s, and t, where the k^ are three independent external momenta, p is the 
loop momentum, are external polarization vectors, and 



are the usual Mandelstam invariants. By dimensional analysis, each numerator term must 
contain four momenta in addition to being linear in all four e^s. We also set ki ■ £j = and 





(3.3) 
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FIG. 5: The two unitarity cuts used to determine the amplitude. The exposed intermediate legs 
are on shell. 

impose momentum conservation with A; 4 = —k\ — k 2 — k 3 and k\ • £4 = — k 2 ■ £4 — k 3 ■ £4. This 
yields 468 terms, each with a coefficient to be determined. 

Our first constraint on the coefficients comes from demanding that the box numerator 
obey the rotation and reflection symmetries of the box diagram. This leaves us with 81 
free coefficients. An ansatz for the full amplitude is then obtained by using the duality 
relations (13. ip . (13 .2p to determine numerators for all other diagrams. 




(4) (5) 




(6) (7) 
FIG. 6: The seven diagrams for the color-ordered amplitude with ordering (1,2,3,4). 

The next step is to determine coefficients in the ansatz by matching to the unitarity cuts 
of the amplitude. It is convenient to use a color-ordered form of the amplitude [48] for this 
matching. The seven diagrams contributing to the color-ordered amplitude, that is the co- 
efficient of the color trace iV c Tr[T ai T a2 T a3 T a4 ], are shown in Fig. O The other color-ordered 
amplitudes are simple relabelings of this one. For the one-loop four-point amplitude, the s- 
and t-channel unitarity cuts — shown in Fig. [5] — are sufficient to determine this color-ordered 
amplitude up to terms that integrate to zero. One straightforward means for determining 
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the cuts is to construct the amplitude in Feynman gauge and then take its unitarity cuts at 
the integrand level prior to integration. This automatically gives us an expression for the 
cuts valid in D dimensions without any spurious denominators (such as light-cone denom- 
inators from physical state projectors). This matching procedure nontrivially rearranges 
the amplitude so that BCJ duality is manifest. After matching the cuts, we also impose 
the fourth auxiliary condition to tame the bubble-on-external-leg contributions. Finally we 
impose that the tadpole numerators vanish. Including all the auxiliary constraints with 
these conditions, we can solve for all but five free coefficients. Because the s- and t-channel 
unitarity cuts are independent of these parameters, the integrated amplitude should not 
depend on them. 

Using the shorthand notation, 

Pi=P, P2=P~h, P3 = P - h - k 2 , p 4 = P + k 4 , 
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and setting the free parameters to zero for simplicity, the box numerator is 



^1234;p 
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Pi 
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7D s -230 
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+ cyclic 



(3.5) 



where D s is a state-counting parameter, so that D s — 2 is the number of gluon states 
circulating in the loop. The notation '+ cyclic' indicates that one should include the three 
additional cyclic permutations of indices, giving a total of four permutations (1,2,3,4), 
(2, 3, 4, 1), (3, 4, 1, 2), (4, 1, 2, 3) of all variables £j, fcj,Pi, s, t. Plain-text, computer-readable 
versions of the full expressions for the numerators, including also gluino and scalar-loop 
contributions, can be found online |49|. In eq. f 1 3 . 5 j) . we have written the expression for 



the box numerator in a different form than available online in order to exhibit the cyclic 
symmetry. 
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FIG. 7: The planar and nonplanar double-box integrals appearing in the identical-helicity pure 
Yang-Mills amplitude. 



We have explicitly checked that after reducing the pure Yang-Mills amplitude to an inte- 
gral basis, 1 the expression is free of arbitrary parameters and in D = 4 matches the known 
expression for the amplitude in ref. j50[, after accounting for the fact that the expression in 
that paper is renormalized. The reduction for four- dimensional external states was carried 
out by expanding the external polarizations in terms of the external momenta plus a dual 
vector [51] . 

As another simple cross check, we have extracted the ultraviolet divergences in D = 6, 8 
and compared them to the known forms. In D = 6, 8, with our fourth auxiliary constraint 
there are no ultraviolet contributions from bubbles on external legs. This allows us to directly 
extract the ultraviolet divergences by introducing a mass regulator and then expanding 
in small external momenta using the methods of ref. j52|. We find complete agreement 
with both earlier evaluations in ref. j53[. We have also compared this to an extraction of 
the ultraviolet divergences directly using dimensional regularization without introducing an 
additional mass regulator and again find agreement. 



B. Two Loops 

As we shall discuss in Section IIVt the four-graviton amplitude in the double-copy theory 
is ultraviolet finite at one loop. To test whether this continues at two loops we need the two- 
loop amplitude. As it turns out, the identical-helicity amplitude is sufficient for our purposes. 
We therefore now turn to finding a form of the two-loop identical-helicity amplitude where 
BCJ duality is manifest. It would be interesting to obtain a general two-loop construction 
valid for all states in D dimensions, but we do not do so here. 

The identical-helicity pure Yang-Mills amplitude has previously been constructed in 
ref. [Hj]. Ref. \^ notes that a representation where the numerators satisfy the BCJ du- 
ality can be obtained by appropriately rearranging the numerators. Here we observe that 
the original representation in ref. [19] also satisfies the duality, up to additional diagrams 
that integrate to zero and are indetectable in ordinary unitarity cuts. 

The two-loop all-plus helicity amplitude for pure Yang-Mills is given by [l9[ 

A[ 2 + , 3 + , 4 + ) = g 6 - ^ [Cr234^i234 + ^12^34^12^34] , (3.6) 



We thank R. Roiban for cross-checking our computation. 
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FIG. 8: The bow-tie integral appearing in the identical-helicity pure Yang-Mills amplitude. 



where the sum runs over all 24 permutations of the external legs. The prefactor 1/4 accounts 
for the overcount due to symmetries. Cf 234 and Cf 2 p M are the color factors obtained from 
the planar and nonplanar double-box diagrams shown in Fig. [7J respectively, by dressing 
each vertex with an f abc and summing over the contracted color indices. A? 234 and ^4f 2 P 34 
are then the planar and nonplanar partial amplitudes. They are given by fl9j 



^1234 



4 NP 
^12:34 



iTTl 



{D s - 2) 



(p + q) 2 X 2 p X 2 q {{p + q) 2 + s) + 16s ((A p • X q ) 2 - X 2 p X 2 q ) 

+ (D s - 2) (5 (a^A 2 + X 2 p X 2 p+q + X 2 q X 2 p+q ) + 4(p + q)\X 2 p + A 2 )(A P • A,))] (s, t) , 
.T,X 4 NP [(D s - 2) (AJAJ + X 2 p Xl +q + AJAJ + J + 16 ((A p • X q f - AjAj)] (s, t) , (3.7) 



where the permutation invariant kinematic prefactor is given by 

r _ [12] [3 4] 
(12>(34>- 

We have rearranged the form of the amplitude given in ref. 1191 very slightly by absorbing the 



(3.8) 



bow-tie contribution depicted in Fig. [S] into the planar double-box contribution in Fig. [7(a) 
by multiplying by a factor of (p + q) 2 in the numerator to cancel the central propagator. 
The planar double-box and nonplanar double-box integrals are 



z*[P{Kp,q,ki)](s,t) 



dp d q 



V(Xi,p,q,ki 



{2tt) d {2tt) d p 2 q 2 {p + q) 2 (p - ktfip - h 



k 2 ) 2 (q - h) 2 (q 



U 2 



dPp d u q 



D, 



V(Xi,p,q, ki 



(2ir) D (2n) D p 2 q 2 (p + q) 2 (p - h) 2 (q - k 2 ) 2 (p + q + k 3 ) 2 (p + q + k 3 + k A ) 



(3.9) 



where X p , X q , and X p+q represent the (— 2e)-dimensional components of loop momenta p, q, 
and (p + q). We note that the planar double-box numerator is nonlocal since it contains a 
term with a factor of 1/s. Nonlocal BCJ representations have been previously found both 
at tree level (3, EH, 12] and loop level 0, |45| . As we now show, the form of the amplitude 
given in eq. (13 .TP manifestly satisfies BCJ duality, once we add an additional diagram which 
vanishes after integration. 

To see if the integrand in eq. (13. 7p satisfies BCJ duality, for each propagator in the planar 
and nonplanar double-box integrals, we consider Jacobi relations (I2.2p around each internal 
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(a) (b) 
FIG. 9: The triangle-box and the bubble-box integrals. 
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FIG. 10: The nontrivial duality relation satisfied by the numerators of the identical-helicity two- 
loop amplitude. The shaded (red) leg marks the central leg of the applied Jacobi identity. 



propagator. Up to relabelings, the only nontrivial duality relation is depicted in Fig. [TUJ 
The other relations shown in Fig. [TT] hold automatically without needing to add additional 
nonvanishing diagrams. In relations (a)-(c), the numerator of the third diagram vanishes, 
while in (d) the invariance of the prefactor T under permutations of external legs and the 
independence of the A's on external momenta make it simple to check that the relation holds. 
Using the diagram numerators from eq. (13. 7ft . the nontrivial duality condition depicted in 
Fig. ITUl implies the existence of a nonvanishing third diagram, shown in Fig. |9^a), which we 
call the triangle-box diagram, 



a tri-box 
/1 2;341 



iT 2l ri_box [4(D a 



+ 



2)(p + q) 2 (X 2 p + X 2 )(X p -X q ) 



(D s - 2f 



(p + q ) 2 \ 2 p \ 2 q (( p + q ) 2 + s )]( s ,t) 



(3.10) 



where 



-7-tri-box 
X 4 



[V(Xi,p,q, ki)](s,t) 
d D p d D q 



V{Xi,p,q, hi 



(2tt) d (2tt) d p 2 q 2 (p + q) 2 (p - ^(p + k 2 + k 3 ) 2 (p + k 2 ) 2 (q - k 2 ) 



(3.11) 



However, this additional triangle-box contribution (I3.10P integrates to zero in dimensional 
regularization and does not contribute to unitarity cuts. To see this, we cancel a factor of 
(p + q) 2 in the numerator with the l/(p + q) 2 propagator. This turns the integral into the 
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FIG. 11: The numerator duality relations for the identical- helicity two-loop amplitude that hold 
automatically. The third diagram in each of (a)-(c) have vanishing numerators and (d) follows 
from simple symmetry considerations of the nonplanar numerators. 



bubble-box integral depicted in Fig. Mjo), 



a tri-box 
/1 2:341 



-jbub-box 



[4( J D S -2)(AJ + A5)(A P -A 
(D 



+ 



2) 2 



X 2 p X 2 q ((p + q) 2 + s)](s,t) 



(3.12) 



where 



d u p d D q 



V(Xi,p,q, h 



(2tt) d {2ti) d p 2 q 2 {p - k^ip + k 2 + k 3 ) 2 (p + k 2 ) 2 (q - h 



\2 ■ 



(3-13) 

The bubble-box integral is a product of two one-loop integrals, one of which is a bubble 
integral with an on-shell massless external leg. Since it is a scale-free integral, it vanishes 
in dimensional regularization, resulting in no contribution to the total amplitude from the 
triangle-box diagram. Moreover, it is not detectable in the standard cuts of the s, t, or 
u = — s — t channels, which is why it is not included in ref. 19 . 
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All other diagrams besides the planar double-box, nonplanar double-box, and triangle-box 
have vanishing numerators by the duality relations. We therefore have a duality-satisfying 
representation for the all-plus helicity pure Yang-Mills amplitude, 



-4i 2) ( 



l + ,2 + ,3+,4 H 



« 6 E 

s 4 



4 < - y 1234 /1 1234 ^ . 



Ly 12;34 /1 



I _/^ytri-box /itri-box 
12:34 "T o°l;234 /1 1;234 



(3.14) 



where the triangle-box contribution integrates to zero. In Section IIVB[ we will use the 
double-copy relation ( 12. 5ft on these numerators to study the ultraviolet behavior of gravity 
coupled to a dilaton and an antisymmetric tensor at two loops. 



IV. ULTRAVIOLET PROPERTIES OF GRAVITY 

We now turn to the ultraviolet properties of the gravity double-copy theory consisting 
of a graviton, dilaton, and antisymmetric tensor, from the perspective of the double-copy 
formalism. The theory generated by taking the double copy of pure Yang-Mills corresponds 
to the low-energy effective Lagrangian of the bosonic part of string theory [HiJ , 

C = ^9 (J^R + + l -e- 2 ^/^H, vp H^ , (4.1) 

where H^ up = d p A up + d u A pp + d p A pu , and A pu = —A u ^ is the rank-two antisymmetric tensor 
field. 



Pure Einstein gravity is one-loop finite in four dimensions [32lj . However, when coupled 
to a scalar (dilaton) [32j or to a rank-two antisymmetric tensor 34 1. the theory is divergent. 
We find that the double-copy theory coupled to both a dilaton and an antisymmetric tensor 
is also divergent, although for all these theories the four-point amplitudes with at least one 
external graviton are finite, as expected from simple counterterm arguments. We will show 
that the cancellation no longer holds at two loops, and the theory has an R 3 counterterm, 
in much the same way as it does for pure Einstein gravity (5o| . In six dimensions, pure 
Einstein gravity is ultraviolet divergent at one loop [35] . We find the same to be true in our 
double-copy theory, and we find a divergence in eight dimensions as well. We will give the 
explicit form of the divergences for these cases. 



A. One Loop 

1. Four Dimensions 

In four dimensions, there is no one-loop four-point divergence when one external leg is a 
graviton l34j because the potential independent counterterms for such divergences vanish 



on shell or can be eliminated by the equations of motion. Using the double-copy formula 
(12. 5p . we have explicitly confirmed finiteness in one- loop four-point amplitudes containing 
at least one external graviton, with the remaining legs either gravitons, dilatons or anti- 
symmetric tensors. We obtain the gravity numerator from the double-copy formula (12. 5p 
by taking the two Yang-Mills numerators, hi and rii, to be equal to the BCJ form of the 
Yang-Mills numerator (13. 5p . As an interesting cross check, we have obtained an asymmetric 
representation of the gravity amplitudes by taking the hi to be the numerators that satisfy 
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BCJ duality and the rij to be numerators obtained by gauge-theory Feynman rules in Feyn 



man gauge, similar to the procedure used recently for half- maximal supergravity [23l. |24| . By 
generalized gauge invariance jl, @| , this should be equivalent to the symmetric construction. 
Indeed, we find identical results for the ultraviolet divergences. 

To evaluate the ultraviolet diver gen ces, we expand in small external momenta to reduce to 



logarithmically divergent integrals 52] . We then simplify tensor integrals composed of loop 
momenta in the numerators by using Lorentz invariance, which implies that the integrals 
must be linear combinations of products of metric tensors rf' v . (See ref. [l3[ for a recent 
discussion of evaluating tensor vacuum integrals.) With the insertion of a massive infrared 
regulator, we finally integrate simple one-loop integrals to find the potential ultraviolet 
divergence. Due to our auxiliary conditions, contributions from bubbles on external legs 
vanish, as they would for ordinary gravity Feynman diagrams. We therefore obtain our 
entire result from box, triangle, and bubble-on-internal-leg diagrams. 

For completeness we have also computed the divergences directly in dimensional regular- 
ization without introducing a mass regulator, using techniques similar to those for two loops 
in Appendix After subtracting the infrared divergence as computed in Appendix [BJ we 
find complete agreement with our result found using vacuum integrals. 

We obtain an expression for the divergence in terms of formal polarization vectors. By 
taking linear combinations of the product of polarization vectors from each copy of Yang- 
Mills, we can project onto the graviton, dilaton, and antisymmetric tensor states. In D = 4 
this is conveniently implemented by using spinor helicity j56[. Graviton polarization ten- 
sors correspond to the 'left' and 'right' copies of Yang-Mills according to — > ^l^Rv 
and — > £^ e^. For the dilaton and antisymmetric tensor, we symmetrize and anti- 
symmetrize in opposite- helicity configurations according to e^ u —> -^{z^^Rv + £ Lfi £ Rv) anc ^ 
£ ^u ~~ >" TJ^l^r.!/ ~~ £ Ln £ Ru)- By substituting the explicit polarizations, we find that all 
configurations where at least a single leg is a graviton are free of ultraviolet divergences, 

M (1) (l\2,3,4) =0, (4.2) 

div. 

where leg 1 is either a positive- or negative-helicity graviton, and the other three states are 
unspecified. 

We however find divergences for the cases with no external gravitons. For the four-dilaton 
amplitude, we find 

*. **>l - 7 (f) 4 (^ n32 'r +3g ' (* 2 + <■ + *> ■ w> 

corresponding to the operator, 

•k\4 1 1132 - 92D, + 3D? 



(4tt) 4 240 



(D^Dy) 2 . (4.4) 



This result is similar to the one obtained long ago by 't Hooft and Veltman [32]. However, 
in our case we have an antisymmetric tensor which can circulate in the loop, altering the 



numerical coefficient. We note that the operator in ref. [32| looks different than above, but 
it can be written in a similar way through use of the field equations of motion. 

The amplitude with four antisymmetric tensors is also one-loop divergent in four di- 
mensions. In four dimensions, the antisymmetric tensor is dual to a scalar field, so we 
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expect the divergence to be the same as that for dilatons. Indeed, the divergence in the 
four-antisymmetric-tensor amplitude for a theory with an antisymmetric tensor coupled to 
gravity is equal to that of the four-dilaton amplitude in a theory of a dilaton coupled to 
gravity j34j. In congruence, we find the divergence for four external antisymmetric tensors 
to also be given by the same expression as the four-dilaton divergence (14.3)) . 

M {1) (1 A ,2 A ,3 A ,A A ) = A** 1 ^!^, 3*4*) • (4-5) 



div. 



div. 



In terms of the antisymmetric tensor fields, the divergence is generated by the operator, 

e \2J (4tt) 4 2160 \ i»>p ) \ ) 

The counterterm that cancels the divergence is given by the negative of this operator. 

In addition to the above divergences, there is also a divergence in the D = 4 
amplitude. This divergence is given by 

1 z /1116 - 76D S - D 2 



_M«(r^,2*,3V A ) 



div. e V27 (4vr) 4 V 120 



s 2 



-1124 + 84L> S -L> 2 /2 ... , 
+ ^ -if +«')), (4-7) 



which corresponds to the operator, 



1 1 /1124 - 84L>, + D 



2 J (4vr) 4 V 60 



1132 - 92£> s + 3D 
360 



H^H^D^D^ . (4.8) 



2. Six Dimensions 



In six dimensions for external gravitons, the only independent invariant operator at one 
loop [57] is 

Rap^R^Rp^. (4.9) 

This corresponds to the known D = 6 one-loop divergence of pure Einstein gravity given 
in ref. j35[. We have computed the coefficient of the D = 6 divergence for the double-copy 
theory of a graviton coupled to a dilaton and an antisymmetric tensor. In this case, the 
divergence is given by the operator, 



1 1 {D s -2 



\2 



RaRnvR^'* 7 Rp^ ■ (4.10) 

e(4vr) 3 30240 ^ p V ; 

Appropriate powers of the coupling are generated by expanding the metric around flat space, 
9iiv = fyiv + ftfyuv Although we do not include the explicit forms of the counterterms here, we 
have also found divergences for the following amplitudes (as well as their permutations and 
parity conjugates) involving external dilatons and antisymmetric tensors, where we restrict 
the external states to four dimensions: 



.M (1) (l ,2 + ,3 + ,4+), _M (1) (l*,2^,3 + ,4+), M {1 \l*, 2*, 3* 4+) , M (1) (l*, 2 , 3* 
.M«(l A ,2 A ,3+,4 + ), M W (1 A ,2 A ,3*,A + ), M {l \l A } 2 A , 3*, 4*) . (4.11) 
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3. Eight Dimensions 



In eight dimensions, there are seven linearly independent R operators [58 



On shell, the combination, 

U = 



T 2 

T 5 
T 7 

T-i 



p pM^P p a p7<5kA 
J^p.vpa^ X "/6 k 



R 



T3HV p^7 



R p*6 



p R/^ po"7<5K 

p P^ P' 5 ^ pcr<57K 
^fiupa-n A7 S k 11 

P pA* P p^ 7 T3V&OK 
-rV/po-n- A 7 5 k 

p pA* P pA v Tj-ySaK 



^ + T 2 - - T 4 + 2T 5 - T 6 + 2T 7 
16 8 



(4.12) 



(4.13) 



is a total derivative, so only six of the Tj are independent on shell. In terms of these 
operators, the divergence for gravity coupled to a dilaton and an antisymmetric tensor at 
one loop in D = 8 is 



1 1 1 
e (4tt) 4 1814400 



[(4274 - 899£> s + HD 2 )^ - 40(466 - 103£> s - 2D 2 S )T 2 
- 2(1886 + 319£> s - D 2 S )T 3 - 180(1034 + D S )T A 
+ 16(1196 + 34£> s - D 2 S )T Q + 64(12454 + 71D S + D 2 )T 7 + cU] 



(4.14) 



where c is a free parameter multiplying the total derivative (14. 13ft . 

We have also found that the following four-point amplitudes involving dilatons and anti- 
symmetric tensors diverge in D = 8: 



/W (1) (l*,2^,3 + ,4 + ) , /Vl (1) (i^2*,3 + ,4-), yVl (1) (l , 2*, 3*, 4*) , M {1) {1 A , 2 A , 3 + , 4+) , 
yW (1) (l A , 2 A ,3 + ,4"), 7 vi( 1 )(i^,2 A ,3^,4 + ), M«(r\ 2 A , 3^) , M«(1 A , 2 A , 3 A , A A ) , 

(4.15) 

where we have again chosen the external states to be four dimensional. The other configu- 
rations are finite. 



B. Ultraviolet Properties of Gravity at Two Loops in Four Dimensions 

Pure Einstein gravity in D = 4 is one-loop finite, but it does diverge at two loops jHHJ. 
This suggests that the two-loop four-graviton amplitude, including also the dilaton and 
antisymmetric tensor, should diverge as well. For external gravitons, the only independent 
operator is the same i? 3 operator for one loop in six dimensions (14. 9p . Our aim is to find its 
coefficient. 
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FIG. 12: The R 3 operator dia grams that contribute to the all-plus helicity four-graviton amplitude. 
The solid dot represents vertices generated by the R 3 operator. 



The R 3 operator generates a nonvanishing four-point amplitude for identical-helicity 
gravitons, illustrated in Fig. [12j This means that we can determine its coefficient by comput- 
ing the four-graviton all-plus helicity amplitude. Fortunately, as we discussed in Section UTTl 
we have the BCJ form of the required all-plus helicity Yang-Mills amplitude. Applying the 
double-copy formula (12.51) to the Yang-Mills amplitude in eq. (13.141) immediately gives us the 
corresponding gravity integrand, simply by squaring the numerators. Thus, for the planar 
double-box, the nonplanar double-box, and the triangle-box, respectively, we have 



<234 



(d s - 2)s (a;aj + a;a; +9 + a^j + ie s ((a, • \f - a^) 



+ A{D s -2){p + q )\Xl + \l){\ p .\ q ) 



+ 



(D 8 



(p + q?>?X((p + <i) 2 + s) 



(s,t) 



j ^-bow-tie 



(D s - 2)s (ApAg + ApAp +g + A^Ap +g ) 

+ 16s ((A p -A,) 2 -A^)) 2 ]( S ,t) 
2 ((D s - 2)s (AjA 2 . + A 2 A 2 + + 16s ((A p • \ q f - AjA 2 ,)) 



j •jbiow-tie 



x ^ 4 ( J D S -2)(A 2 + A 2 )(A P -A 9 ) + 
(p + qf{A(D s -2){\l + \l){\ p .\ q ) 



{Ds - 2) 



-2 „2 T NP 



(D s - 2) (AjA 2 . + AjAj +g + AJA; + J 

+ 16((A p -A g ) 2 -A^)) 2 l( S ,t) 
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]\ yftvi— box 
JW 2;341 



ri— box 



2){ P + qf{\l + \ 2 q ){\ p -\ q ) 



+ 



2 -7-bubb— box 



, {D s - 



[s,t) 



KK((P + qV + s] 



2 ~i 



[s,t). 



(4.16) 



We have explicitly confirmed that s-, t-, and w-channel unitarity cuts are satisfied. We did 
so numerically keeping the internal states in integer dimensions D = 6 and D = 8. 

To obtain the ultraviolet divergences, we integrate the amplitudes in dimensional regu- 
larization. As was the case in pure Yang-Mills, the triangle-box does not contribute. It once 
again reduces to a bubble-box integral, which vanishes in dimensional regularization. We 
carry out the extraction of the ultraviolet divergences in two ways. In the first approach 
we simply use dimensional regularization and then subtract the known infrared divergences, 
leaving only the ultraviolet ones. In the second approach we introduce a mass regulator 
to separate the ultraviolet singularities from the infrared divergences, as carried out in Ap- 
pendix Either method yields the same result. In fact, the second method also shows that 
the bubble-box integral is not ultraviolet divergent. Either way, we therefore have that the 
bubble-box does not contribute to the divergent parts, 



A^ 2 )(l + ,2+,3 + ,4^ 



div. 



1 



Si 



div. 



(4.17) 



The dimensionally regularized integrals are performed in Appendix [A] Eq. flAlOO gives 
the planar double-box integrals; eq. (1A17j) gives the nonplanar double-box integrals; and 
eq. (1A20P gives the bow-tie integrals. The infrared divergence from Appendix [B] is 



A4 (2) (l + ,2+ 3+4+) 



IR div. 



1 fK\ 6 

~e\2) (47r^ 4 



'—T 2 (jPs 2)2 (s 2 + t 2 + u 2 ) 

^4 12 Q V J 



X 



S log(-^j+tlog ^j+„log^ 



(4.18) 



We insert the divergent parts of the integrals evaluated using dimensional regularization 
into eq. (I4.16p . then insert these results into eq. ( I4.17P and perform the permutation sum. 
Finally we subtract the infrared divergence and arrive at the two-loop ultraviolet divergence 
of gravity coupled to a dilaton and an antisymmetric tensor for four external positive-helicity 
gravitons: 



A^ 2 )(l+,2+,3 + ,4^ 



lUV div. 



k\ 6 i 
2) (4tt) 4 



T 2 (47£ 2 



748ZX + 2268Ww 



240 



(4.19) 



For our second method, we evaluate the ultraviolet divergences of the required integrals by 
going to vacuum and using a massive infrared regulator, sidestepping the need to subtract the 
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infrared divergence. The ultraviolet divergences of the individual integrals are calculated in 
Appendix O After permutations, the contributions of the planar and nonplanar components 
are 

jw p (i+,2+,3+,4+)L r ,. - 1 (1)\,\ A r 2 {l1Dl " lsgDs + m)stu 



uvdiv. e \ 2 J (4tt) 4 60 

v ' ' ' ;iuv div. e V2/ (4tt) 4 240 v ; 

Summing these two contributions, we find complete agreement with eq. (I4.19p . 

We can re-express the two-loop divergence in terms of the operator that generates it. By 
matching the amplitude generated by the diagrams with an R 3 vertex shown in Fig. [12] to 
the divergence in eq. (14.191) . we find that the operator, 

1 (K _y 1 47^- 748/^ + 2268 , 

e\2J (4vr) 4 14400 ^ V ; 

generates the two-loop divergence for gravity coupled to a dilaton and an antisymmetric 
tensor. 



V. CONCLUSION 

In this paper we constructed a representation of the one-loop four-point amplitude of 
pure Yang-Mills theory explicitly exhibiting the duality between color and kinematics. This 
construction is the first nonsupersymmetric example at loop level valid in any dimension 
with no restriction on the external states. The cost of this generality is relatively complicated 
expressions in terms of formal polarization vectors. 

The duality between color and kinematics and its associated gravity double-copy structure 
has proven useful for unraveling ultraviolet properties in various dimensions [E, [l3|, 23. 24, 31 



Using the one-loop four-point pure Yang- Mills amplitude with the duality manifest, we ob- 
tained the integrand for the corresponding amplitude in a theory of a graviton, dilaton, and 
antisymmetric tensor. In D = 4, we found that one-loop four-point amplitudes with one or 
more external gravitons are ultraviolet finite, while amplitudes involving only external dila- 
tons or antisymmetric tensor fields diverge. This result is similar to those of earlier studies 
involving gravity coupled either to a scalar, an antisymmetric tensor, or other matter and 



is in line with simple counterterm arguments |32M34j . We gave the explicit form, including 



numerical coefficients, for all four-point divergences in this theory. Since our construction 
is valid in any dimension, we also investigated the ultraviolet properties of the double-copy 
theory in higher dimensions. In particular, we showed that in D = 6, 8 the one-loop four- 
graviton amplitudes diverge, as expected, and gave the explicit form of these divergences 
including their numerical coefficients. 

In order to investigate whether the observed D = 4 ultraviolet finiteness of the amplitudes 
with one or more external gravitons continues beyond one loop, we also computed the coeffi- 
cient of the potential two-loop R 3 divergences. This is greatly simplified by the observation 
that an R 3 divergence can be detected in the identical-helicity four-graviton configuration. 
The required gravity amplitude was then easily constructed via the double-copy property, 
by first finding a representation of the pure Yang-Mills amplitude that satisfies the duality. 
The existence of such a representation has already been noted in ref. [|3j]. Here we used a 
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slightly different one, matching the original form of the two-loop identical-helicity ampli- 
tude given in ref. [l9j], but adding an additional diagram that integrates to zero. We found 
that the two-loop amplitude with external gravitons is indeed divergent and that the R 3 
counterterm has nonzero coefficient. This is not surprising given that pure Einstein gravity 
diverges at two loops [55j ]. Our paper definitively shows that, as one might have expected, 
the double-copy property by itself cannot render a gravity theory ultraviolet finite. For 
ultraviolet finiteness, an additional mechanism such as supersymmetry is needed. Further 
progress in clarifying the ultraviolet structure of gravity theories will undoubtedly rely on 
new multiloop calculations to guide theoretical developments. We expect that the duality 
between color and kinematics will continue to play an important role in this. 
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Appendix A: Dimensionally Regularized Integrals 



In this appendix, we explicitly compute the divergent parts of dimensionally regularized 
two-loop integrals in D = 4 — 2e, appearing in Section TlV Bl In general, both ultraviolet and 
infrared divergences appear as poles in e so we must subtract the infrared ones in order to 
obtain the ultraviolet ones. 

We start with the planar double-box integral, displayed in Fig. [7(a), following the dis- 



cussion in Section 4 of ref. [59 

li[P{\ h p,q, ki)](s,t) 
d D p d D q 



V(Xi,p,q,ki 



{2tt) d {27r) D p 2 q 2 {p + q) 2 {p - kxfip -k x - k 2 ) 2 {q - h) 2 {q - k 3 - fc 4 ) 2 ' 



(Al) 



Using Schwinger parameters, we rewrite the planar-double box integral with constant nu- 
merator as 



2? [!](«,*) 



(4vr 



-ri 

w\D 11 



1=1 



dU [A p (T)]-t exp 



Qp(s,t,U 
A P (T) 



where 
and 



A P (T) = (T p T q + TpTpg + TgTj 



pqj 5 



T 

J. . „ 



f 3 + U + t 



5 ■ 



T 



U+U + tn 



T 



/"/ 



(A2) 



(A3) 



(A4) 



T p , T q , and T pq are sums of Schwinger parameters corresponding to propagators with loop 
momenta p, q, and p + q respectively. We also have 



(A5) 



Q P (s,t,ti) = -s tit 2 T p + t 3 UT q + t 6 (ti + t 3 )(t 2 + U) -tt 5 t 6 t 7 
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To account for factors of A p , A 2 ,, and X 2 +q in the numerator, we take derivatives on the 
(— 2e)-dimensional part of the (Wick- rotated) integral: 

J dX p 2 ^dX q 2 %xp [-T P X 2 - T q \ 2 q - T m X 2 p+q ] oc [A P (T)] £ , (A6) 

with respect to T p , T q , and T pq . This introduces additional factors to be inserted in the 
integrand in eq. (1A2I) . For example, 

(AJ) 4 -6(l-6)(2-6)(3-6) 

(A 2 ,)^ 2 . -+ e 2 (l - e)(2 - e ) (T ^ p + ( ^ 3 )2 - 6(1 - e)(2 - e)(3 - e) ^^4^ , 
(A^A^A 2 e 2 (l - e) 2 -^— + C (l - e)(2 - e) fS±^±^ (A7) 



-e(l- e )(2- e )(3-e)— Vti- 
V A a ; A P (T) 4 

We account for extra factors of A P (T) by shifting the dimension D — )■ D — 2a. Following 
Smirnov [60j] , we change six of the seven Schwinger parameters to Feynman parameters with 
the delta-function constraint Ylii^& a i = 1 ; 



X 4 p [P(A p ,A g )]( g ,t) = r[7 r . ^ + 7] f'daJl ^ {T }\^ +1 D{a 

(A8) 

where D(a.j) represents the extra factors in one term of eq. fl ATI) . with t% — >■ a«. The 
parameter 7 counts the factors of in .D and can take on values 0, 2, and 4 for the 
integrals under consideration here. Next we perform a change of variables that imposes the 



1/ 1 



delta-function constraint 60 



ai = #i.6> «2 = (1 - 6>)(1 - &) , o 3 = /3 2 £i, 04 = ^5(1-^2), 

«5 = /3 2 (l-6), a 7 = A(l-6), & = & = £b6- (A9) 

We then integrate these parameters to obtain a Mellin-Barnes representation, which we 
again integrate. Finally we arrive at the dimensionally regularized results of our required 
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planar double-box integrals: 



Xf[(A 2 ) 4 ]M 



Xf[(Aj +q ) 4 ](M 



Xf[(Aj) 3 AjM 
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(All) 



All integrals above are symmetric under A p -H- A g . 

Next we look at the nonplanar double-box integrals: 



l^[V{\ hP ,q,h)]{s,t) 
d D p dPq 



V(\i,p,q, h 



(2n) D (2n) D p 2 q 2 (p + q) 2 (p - h) 2 (q - k 2 ) 2 (p + q + k 3 ) 2 (p + q + k 3 + k A ) 2 ' 

(A12) 

whose evaluation follows that of the planar double-box integrals quite closely. Anp(T) takes 
the same form as A P (T) in eq. (1A3[) . except that 

T p = ti + t 2 , T q = t 3 + , T pg = t 5 + t 6 + t 7 . (A13) 

We then also have 

QnpO, t, u, U) = S (t!t 3 t 5 + t 2 Ut 7 + t 5 t 7 (T p + T q )) - t t 2 t 3 t 6 - u txUU . (A14) 
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In this case, we find it advantageous to only change the four Schwinger parameters associated 
with T p and T q to Feynman parameters, resulting in 

zrwvA,)] = t[7 ^d i] ( Al5 ) 



X 



n/ dcull detail [ NP )l 7 _ g+7 D(aO. 
7=5 ^ 7=1 -J" \ 1=7/ [Qnp(s, 



We impose the delta-function constraint via further redefinition: 

ai = 6(i -6), « 2 = 66, a 3 = (1-&X1-6), a4 = (1-6)6- (A16) 

Once again we can straightforwardly integrate the parameters and use the Mellin-Barnes 
representation to evaluate our required nonplanar double-box integrals: 

X^\(X^](st)-X' NP 1 215fi2 ± M2st ± M2t2 | C( C Q) 
X 4 NP [(Aj + /]( S ,t) = -^ I -^ + O( e ), 



p+qJ n ' ; (4vr) 4 80e 
X 4 NP [(Aj) 3 Aj + J( S ,t) = G( e °), 



4 lv P ; v q j jv , ; ^ 4 25200se v ; ' 



X 4 NP [(A;) 2 (AJ + /](M) = ^ I ^ + O(e ), 

xr[(Aj) 2 A^ ( M) = ^ I ^ + 0( e °), 
xr[A5A;(A^) 2 ]( S ,t) = 0( e °), 

xr P [A5(A^) 3 ]( S ,t) = + O(e ) , (A17) 
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+ O(e ) . (A18) 



As with the planar results, the above are valid under the exchange \ p -H- X q . 
Finally we evaluate the bow-tie integrals: 

d D p d D q V{Kp,qM) 



I 



(2w) D (2%) D p 2 q 2 (p - h) 2 (p — ki — k 2 ) 2 (q - k 4 ) 2 (q - k 3 - h) 2 ' 



(A19) 
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The bow-tie integrals are relatively simple because they are products of two one-loop inte- 
grals. Similar techniques involving Schwinger parameters and Mellin-Barnes representations 
can be used on each one-loop integral. Since bubbles with a massless leg vanish in dimen- 
sional regularization, the replacement (p + q) 2 — > 2p ■ q is valid in the numerator. We also 
use the tensor reduction (X p ■ \ q ) 2 — > A 2 A 2 /(— 2e). The relevant bow-tie integrals are then 
given by 

X4 bow-tic [(A 2 )2(Ap . A(?)2](s) 
J bow-tie [(A 2 )2(Ap . Ag)2(p + g)2](S)t) 

^° w " tic [A^(A p -A 9 ) 2 (p + g) 2 ]( S ,t) 
These are also symmetric under the exchange A p -H- X q . 



1 s 2 

' + O(e ) 



(4tt) 4 720e 
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1 s 2 (10s-t) 
(4tt) 4 15120e 

1 s 2 (12s-t) 
(4tt) 4 28800e 



+ O(e ) 



+ 0(e°) . (A20) 



Appendix B: Two-Loop Infrared Divergence 



In this appendix we obtain the two-loop infrared divergence for the four-point all-plus 
helicity graviton amplitude in the theory of gravity coupled to a dilaton and an antisym- 
metric tensor using dimensional regularization in D = 4 — 2e. We subtract these infrared 
divergences from the total divergence to obtain the ultraviolet divergence. Infrared diver- 
gences in gravity can be obtained by exponentiating the divergence found at the one-loop 



order [16l. l6ll. |62|. In the cases where there is a divergence at one loop, the infrared singular- 
ities are 'one-loop exact'; however, in the all-plus helicity gravitons case, the first divergence 
occurs at two loops. Nevertheless, the same principles apply. More specifically we are con- 
cerned with the exponentiation of the gravitational soft function, which describes the effects 
of soft graviton exchange between external particles. 

Following the discussion of ref. [62|, a gravity scattering amplitude can be written as 



(Bl) 



where S n is the infrared-divergent soft function and H n is the infrared-finite hard function. 
Each quantity in eq. fIBlj) can be written as a loop expansion in powers of (fi;/2) 2 (47re _7B ) € : 



S n — 1 + H n — Hjf\ 



L=0 L=l L=0 

The soft function is given by the exponential of the lowest order infrared divergence: 

2 1 
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An L-loop amplitude can then be written as 
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1=0 v 
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(Jr. 
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(B4) 



For four point amplitudes, we have 
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and the one-loop infrared divergence is given by 
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(B6) 



We used this to subtract the infrared divergence from our dimensionally regularized one- 
loop result in Section IIV A II to isolate the ultraviolet divergence. The four-point two-loop 
infrared divergence is given by 



MP 
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-0-4" 


IR div. 2 


. e . 



(0) , ..(1) 



Mr + —m 



(B7) 



For the all-plus helicity gravitons case, the tree amplitude M^ vanishes. The one-loop 
amplitude can be computed using the double-copy procedure in Section IIV A II and is given 
by H 



A4 (1) (l + ,2+,3 + ,4+ 



.2) (4vr) 2 V(12)(34) J 
The two-loop infrared divergence is then 
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Appendix C: Ultraviolet Divergences from Vacuum Integrals 

In this appendix we compute the ultraviolet divergences of the integrals in Section IIVBI 
The techniques are very similar to those used to study the one-loop ultraviolet proper- 
ties of gravity in Section IIVI However, before we can use them, we must deal with the 
-2e)-dimensional components \„, \ q , and X p+q in the numerators of the integrals using the 



techniques in Section 4.1 of ref. [59 



The effect of inserting factors of X p , X q , and \ p + q into the planar and nonplanar double-box 
integrals is very similar to inserting factors of v ■ p, v ■ q, and v ■ (p + q), where 
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Example parameter insertions for factors of A j are given in eq. (1A7I) . For polynomials in t> • p 
and v ■ q, we have 



A J 



(f -p) 6 (f ■ q) 2 



-pq) 
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(C2) 



These are valid for both the planar and nonplanar double-boxes provided the corresponding 
definitions for A, T p , T q , and T pq given in Appendix |A] are used. 

We can also relate polynomials in v -p and v-q to the Aj. The four- dimensional component 
of the loop momenta p can be written as 



(C3) 



where 
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We therefore have 



p 2 + A 2 = p [4 ] • P[ 4 ] = sc?c^ + £c|c| + uc?cg - -stw(c^) 2 



(C5) 



or 



■ — (vp) 2 + V p 
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where 



V p = — p 2 + sd[c^ + tc^Cg + wc^cf . 



(C7) 



Similarly, we have 
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AL, = ~(v ■ (p + q)) 2 + P pq , 



stu 



(C8) 
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where 



V q = — q 2 + s c\c q 2 + t c\c\ + u c\c\ , 

V pq = - (p + ?) 2 + s(c? + c?)(4 + c!) + t(ci + c2)(cg + cl) + «(c? + c?)(cf + 4) . (C9) 

These relations, along with the parameter replacements eqs. ( 1 A7I) . (1C2|) . allow us to rewrite 
the integrals involving factors Aj in terms of integrals involving tensor products between the 
loop momenta and the external momenta. For a general function f(p ■ ki,q ■ hi), we have 
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(CIO) 



where a factor 1/A indicates that a shift in dimension of the integral should be made: 
D— > D + 2, e— > e — 1 (e's in prefactors in eq. ( IClOj) should not be shifted, however). 

Once we have integrals in a form involving tensor products between the loop momenta 
and external momenta, we expand in small external momenta to reduce to logarithmically 
divergent integrals, just as we did in the one-loop case. This gives us vacuum integrals. 
We then reduce the tensors involving loop momenta using Lorentz covariance and insert an 
infrared mass regulator. By integrating we obtain the ultraviolet divergences. Since every 
prefactor in eq. (1C10I) contains a factor of e, to get the ultraviolet divergence, we only need 
the 1/e 2 pole of the integrals on the right-hand side. These leading contributions have no 
dependence on the mass regulator, so we are unaffected by subdivergence issues due to the 
mass regulator. The ultraviolet divergences of the planar and nonplanar double-box integrals 
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The bow-tie integrals do not contain infrared divergences, and their ultraviolet divergences 
were computed in Section |A] Combining all the pieces then gives us the ultraviolet divergence 
in eq. 
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